Anomalous magneto-oscillations and spin precession 
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A semiclassical analysis based on concepts developed in quantum chaos reveals that anomalous 
magneto-oscillations in quasi two-dimensional systems with spin-orbit interaction reflect the non- 
adiabatic spin precession of a classical spin vector along the cyclotron orbits. 



If in a solid the spatial inversion symmetry is broken, 
spin-orbit (SO) interaction gives rise to a finite spin split- 
ting of the energy bands even at magnetic field _B = 0. 
In quasi two-dimensional (2D) systems this B ~ spin 
splitting is frequently analyzed by measuring the magne- 
toresistance oscillations at small magnetic fields i? > 0, 
known as Shubnikov-de Haas (SdH) oscillations. Follow- 
ing a semiclassical argument due to Onsager ||] it has 
long been assumed ||, ||], that the frequencies of 
these oscillations are proportional to the unequal occu- 
pations N± of the spin-split subbands, 

7V± = (e/2^?i)/fH, (1) 

where e is the electron charge and Ti is Planck's constant. 
Recently, experiments and numerical quantum mechan- 
ical calculations have shown that, in general, these os- 
cillations are not simply related to the B — Q spin- 
subband densities However, it has remained un- 

clear when and why Onsager's semiclassical argument 
fails. Here we use a semiclassical trace formula for parti- 
cles with spin, which was only lately developed |^] in 
the context of quantum chaos, in order to show that the 
anomalous magneto-oscillations reflect the nonadiabatic 
spin precession along the cyclotron orbits. Currently 
great efforts are made to obtain a deeper understand- 
ing of spin-related phenomena in semiconductor quan- 
tum structures, in particular due to possible applications 
in spintronics Q. While spin is a purely quantum me- 
chanical property with no immediate analogue in clas- 
sical physics, the present analysis reveals that our un- 
derstanding of spin phenomena can be greatly improved 
by investigating equations of motion for a classical spin 
vector. 

In the presence of a magnetic field perpendicular to the 
plane of a 2D electron system the electrons condense in 
highly degenerate Landau levels that are regularly spaced 
in energy. With increasing field B these Landau levels 
are pushed through the Fermi surface causing magneto- 
oscillations which reflect the oscillating density of states 
(DOS) at the Fermi energy Ep, see, e.g., Q. Onsager's 
semiclassical analysis of magneto-oscillations was based 
on a Bohr-Sommerfeld quantization of cyclotron orbits. 
However, for systems with spin there is no straight- 



forward generalization of Bohr-Sommerfeld quantization 
. The Gutzwiller trace formula provides an al- 
ternative and particularly transparent semiclassical inter- 
pretation of magneto-oscillations that is applicable even 
in the presence of SO interaction. Rather than giving in- 
dividual quantum energies, the trace formula relates the 
DOS of the quantum mechanical system to a sum over 
all periodic orbits of the corresponding classical system. 
As a function of energy E the individual terms oscillate 
proportional to cos[S{E)/h] where S{E) is the action of 
the orbit. 

First we briefly discuss magneto-oscillations of elec- 
trons with effective mass to* in a 2D system without 
SO interaction. Here, the sum over fc-fold repetitions 
of the classical periodic cyclotron orbits corresponds to 
a Fourier decomposition of the DOS as a function of 
the energy E, where the action of the /c-fold revolu- 
tion corresponds to the /cth harmonic 2nkE / (tku) — 
2Trkm* E / (heB) of the DOS with cyclotron frequency 
uj — eB/m*. Thus we see that the DOS for a fixed 
energy E — Ep oscillates as a function of the recipro- 
cal magnetic field 1/i? which is the origin of magneto- 
oscillations. In particular, we get Onsager's formula from 
the first harmonic k — 1. Longer orbits fc > 1, giving rise 
to higher harmonics in the oscillating DOS, are exponen- 
tially damped for small but nonzero temperatures [ p^ so 
that here it suffices to consider A: = 1. In general, the 
Gutzwiller trace formula is an asymptotic relation that 
holds in the semiclassical limit fi ^ Q. However, in the 
particular case discussed above it is an identity. 

Recently it has been shown ^ that in leading semi- 
classical order the SO interaction results in weight factors 
2cos(A:a/2) for the orbits in the trace formula where the 
angle a characterizes the spin precession s = s x 3 oi a 
classical spin vector s along the classical orbit ||l^ . Here 
s is the time derivative of s and 3 is an effective magnetic 
field including the contributions of both SO coupling and 
the Zeeman interaction due to the external magnetic field 
B felt by the spin along the orbit. After k periods of the 
cyclotron motion the spin vector s has been rotated by 
the angle ka about an axis n, see Fig. |l|. We remark that 
the axis n (but not a) depends on the starting point of 
the cyclotron orbit. Like the effective field 3, the angle a 
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FIG. 1: (color) Classical spin precession (bold green line) 
about the effective field B (bold red line) along a cyclotron 
orbit (black) for a GaAs QW. Thin lines represent the mo- 
mentary vectors of the effective field B (red) and the spin s 
(green) along the cyclotron orbit. The momentary vectors for 
3 are normalized with respect to the maximum of |B| along 
the orbit. In the starting point we have chosen s \\ 3. After 
one cycle the motion of the spin vector can be identified with 
a rotation by an angle a about an axis n, as shown in the 
blow-up on the left. Initial and final directions of the spin 
vector s are marked in blue. The system is a 100 A wide 
GaAs-Alo.sGao.s As QW grown in the crystallographic direc- 
tion [113] with 2D density A'^ = 5 x 10^^ cm~^ in the presence 
of an electric field -Ex = 100 kV/cm and a magnetic field 
B = 0.05 T. 



depends on the external field B. It contains a dynamical 
as well as a geometric phase similar to Berry's phase [ p^ . 
Apart from higher harmonics the oscillating part of the 
DOS at the Fermi energy Ep is proportional to 



cos(a/2) cos[2T:m*EF/{heB)] 



(2) 



We have analyzed magneto-oscillations for quasi 2D 
electron systems in semiconductors such as GaAs where 
we have two contributions to the SO coupling. The Dres- 
selhaus term reflects the bulk inversion asymmetry of 
the zinc blende structure of GaAs. If the inversion sym- 
metry of the confining potential of the quasi 2D system 
is broken, we get an additional SO coupling given by the 
Rashba term [|l^. While the Dresselhaus term is fixed, 
the Rashba SO coupling can be tuned by applying an 
electric field Ej^ perpendicular to the plane of the quasi 
2D system [pi . 

In Fig. 0^ we compare the Fourier spectra of the 
magneto-oscillations of the DOS calculated by means of 
a diagonalization of the quantum mechanical Hamilto- 
nian with the spectra obtained from Eq. (^) based on 
an integration of the classical equations of motion for 
the precessing spin. We consider here a 2D electron sys- 
tem in a 100 A wide GaAs- Alp. sGao.sAs quantum well 
(QW) grown in the crystallographic direction [113] with 
constant total density = A^^^. -h iV_ = 5 x 10^^ cm"^ 
and with varying For comparison, the circles mark 
the peak positions which one would expect according to 
Eq. (|l|) for the spin subband densities N± calculated 
quantum mechanically at i? = 0. The Fourier spectra 
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FIG. 2: (a) Quantum mechanical and (b) semiclassical Fourier 
spectra for different values of the electric field for a 2D 
electron system in a 100 A wide GaAs-Alo.sGao.s As QW 
grown in the crystallographic direction [113] with constant 
total density N = 5 x 10^^ cm~^. The open circles show the 
expected Fourier transform peak positions {2-kTi/ e)N± accord- 
ing to the calculated spin subband densities N± at B = 0. 



are in strikingly good agreement. On the other hand, the 
peak positions deviate substantially from the positions 
expected according to the B = Q spin splitting. In par- 
ticular, the semiclassical analysis based on Eq. (|^) repro- 
duces the central peak that is not predicted by Eq. (|^). 
The asymmetry in Fig. ^ with respect to positive and 
negative values of reflects the low-symmetry growth 
direction [113] (Rcf. g). 

An analysis of the classical spin precession along the 
cyclotron orbit reveals the origin of anomalous magneto- 
oscillations. The spin-split states at i? = correspond to 
flxing the direction of spin parallel and antiparallel to the 
effective fleld 3{p) along the cyclotron orbit, where p is 
the kinetic momentum. However, in general the precess- 
ing spin cannot adiabatically follow the momentary fleld 
B{p). This can be seen in Fig. |] where we have plotted 
the momentary field B{p) as well as the precessing spin 
s along a cyclotron orbit. Both the direction and the 
magnitude of 13 change along the orbit. In particular, 
the Dresselhaus term reverses the direction of B when 
\B\ has a minimum. A spin vector that is no longer par- 
allel or antiparallel to B implies that the system is in a 
superposition of states from both spin subbands so that 
the magneto-oscillations are not directly related to the 
-6 = spin splitting. 

For the spin, in order to be able to follow the mo- 
mentary field B{p) adiabatically, the orbital motion must 
be slow compared to the motion of the precessing spin, 
i.e., we must have B <^ \B{p)\ for all points p along 
the cyclotron orbit. Therefore, it is the smallest value 
Bmin — min |iB(p)| along the cyclotron orbit which deter- 
mines whether or not the spin evolves adiabatically. This 
is illustrated in Fig. |l| where the parameters were chosen 
such that initially the spin is parallel to the effective field 
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13 . First s can follow B, but after a quarter period of the 
cyclotron orbit the effective field B reaches its minimum 
Bmin and s starts to "escape" from B. Subsequently, 
the spin vector s is no longer parallel to B also in those 
regions where B becomes large again. We remark that 
adiabatic spin precession does not imply a — but only 
that the rotation axis n is approximately parallel to the 
initial (and final) direction of the effective field B. 

For many years, anomalous magneto-oscillations have 
been explained by means of magnetic breakdown ^ . 
Underlying this approach is a rather different semiclas- 
sical picture where each spin-split subband is associated 
with an energy surface with separate classical dynamics. 
In our treatment, on the other hand, there is only one en- 
ergy surface complemented by the dynamics of a classical 
spin vector. It is the essential idea within the concept of 
magnetic breakdown that in a sufficiently strong external 
magnetic field B electrons can tunnel from a cyclotron 
orbit on the energy surface of one band to an orbit on the 
energy surface of a neighboring band separated from the 
first one by a small energy gap. For spin-split bands the 
separation of these bands is proportional to the effective 
field B, i.e., magnetic breakdown occurs most likely in 
regions of a small effective field B. This approach im- 
plies that the anomalous magneto-oscillations are essen- 
tially determined by the breakdown regions only. (These 
breakdown regions can be identified with mode conver- 
sion points [^.) We want to emphasize that here our ap- 
proach differs fundamentally from these earlier models: 
In the present ansatz spin continuously precesses along 
the cyclotron orbit, i.e., the angle a in Eq. (^ is affected 
by the nonadiabatic motion of s in the regions of both 
small and large B (see Fig. |^) . 

In the adiabatic regime B <C B^in the angle a is given 
by the integral of the modulus of the momentary field B 
along the orbit plus a Berry phase jl^ , 

a = / \B\dt + aB, (3) 
Jo 

where T = 27r/a; is the period of the cyclotron motion. In 
the limit of small external fields {B 0) the Berry phase 
as converges towards a constant and the integrand in 
Eq. (H) can be expanded with respect to a small Zeeman 
term, \B\ « Bo + BiB, where the coefficients Bq and Bi 
are T-periodic in time. Thus in the limit of small external 
fields we obtain a{B) w ao/B+ai with constants ao and 
ai independent of i?. Inserting the last relation in Eq. (||) 
we thus retrieve Eq. (||), i.e., only in the fimit of adiabatic 
spin precession magneto-oscillations are directly related 
to the B = spin splitting. By changing the crystal- 
lographic growth direction of the QW, it is possible to 
tune the value of i3min- In particular, for a QW grown in 
the crystallographic direction [110] the Dresselhaus term 
vanishes for p parallel to the in-plane directions [001] 
and [OOT]. Thus for a symmetric QW without Rashba 
SO coupling we have |B(p)| = B for these values of p 



which implies that there is no adiabatic regime and one 
always observes anomalous magneto-oscillations. 

For a system with Rashba SO coupling but no Dressel- 
haus term both the classicsal equations of motion for the 
precessing spin and the quantum mechanical problem can 
be solved analytically. Here the Gutzwiller trace formula 
exactly reproduces the quantum mechanical density of 
states for i? > 0. For this system, the effective magnetic 
field B{p) along the cyclotron orbit has a constant mag- 
nitude B and for small external fields B ^ the exact 
solution turns into the adiabatic solution so that we get 
no anomalous magneto-oscillations in agreement with an 
earlier quantum mechanical analysis [^. 

Finally, we note that the concepts developed here are 
rather general and, in particular, are not restricted to 
spin- 1/2 systems. Indeed, an analogous semiclassical 
analysis can be carried out for any system with (nearly) 
degenerate subbands. These bands can be identified with 
a single band with a SO coupling acting on an effec- 
tive spin degree of freedom similar to Lipari and Bal- 
dareschi's treatment [|l^ of the multiply degenerate va- 
lence band edge in semiconductors with diamond or zinc 
blende structure. In particular, we expect that our ap- 
proach can be applied to the interpretation of de Haas- 
van Alphen experiments on ultrahigh-purity samples ]2C|] 
that had called in question the established concepts of 
magnetic breakdown. 
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